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Abstract 


A multi-objective optima] design methodology is developed for light- weight, low-cost 
composite structures of improved dynamic performance. The design objectives include 
minimization of resonance amplitudes (or maximization of modal damping), weight, and 
material cost. The design vector includes micromechanics, laminate, and structural shape 
parameters. Performance constraints are imposed on static displacements, dynamic am- 
plitudes, and natural frequencies. The effects of damping on the dynamics of composite 
structures are incorporated. Preliminary applications on a cantilever composite beam illus- 
trated that only the proposed multi-objective optimization, as opposed to single objective 
functions, simultaneously improved all objectives. The significance of composite damping 
in the design of advanced composite structures was also demonstrated, indicating that de- 
sign methods based on undamped dynamics may fail to improve the dynamic performance 
near resonances. 

Nomenclature 


A 

[C),[c] 

E 

f 

u 


Area. 

Global and modal damping matrices. 
Normal modulus. 

Frequency. 

Damped frequency. 
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F(z) 

Objective functions. 

G 

Shear modulus. 

G(z)f 

Inequality constraints. 

h 

Thickness. 

H(z) 

Equality constraints. 

mm 

Global and modal stiffness matrices. 

k 

Volume ratio. 

[M], [m\ 

Global and modal mass matrices. 

P,P 

Global and modal excitation force. 

q 

Modal vector. 

t 

Time. 

U 

Displacement vector. 

u 

Dynamic amplitude. 

V 

Weighting coefficients. 

Vf 

Average fiber volume per unit area. 

tu, 

Stored and dissipated specific strain energies 


Stored and dissipated strain energies. 

z 

Design vector. 

6 

Fiber orientation angle. 

V 

Poisson’s ratio. 

P 

Mass density. 

a 

Stress. 

ip 

Specific damping capacity. 

[*] 

Modal matrix. 

Subscripts 


e 

Finite element. 

f 

Fiber. 

l 

Ply (on-axis). 

L 

Laminate. 


tBold characters indicate vectors. 
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771 

Matrix. 

n 

n-th mode. 

d 

Dynamic. 

Superscripts 

T 

Resonance, 

L 

Lower. 

U 

Upper. 

Direction 

11 

Normal longitudinal. 

22 

Normal in-plane transverse. 

12 

Shear in-plane. 

n 

Normal (isotropic material). 

s 

Shear (isotropic material). 

x,V,z 

Structural axes. 


Introduction 

Fiber composite materials are widely used in structural applications requiring high 
stiffness- to- weight and strength- to- weight ratios, as they readily provide high specific mod- # 
uli, high specific strengths, and tailorable anisotropic elastic properties. A spectrum of 
design methods has been developed, ranging from basic to optimal tailoring of compos- 
ite laminates and structures, in order to meet requirements for high performance and 
light-weight. However, most such design methodologies are primarily based on stiffness 
and strength tailoring, hence, they can predict good static structural performance, but 
not necessarily an optimal dynamic performance as they neglect an important dynamic 
parameter, the damping of composite materials. 

Polymer matrix composites are known to exhibit significantly higher damping com- 
pared to most common metals. The previously stated requirements for advanced light- 
weight structures virtually exclude most traditional sources of passive damping, therefore, 
the option to utilize the damping capacity of polymer-matrix composites appears very at- 
tractive. Reported research on the damping of unidirectional composites and laminates 1-6 
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has shown that the damping of composites is highly- tailorable and is primarily controlled 
by constituent parameters (fiber/matrix properties, fiber volume ratio), and laminate pa- 
rameters (ply angles/thicknesses, stacking sequence). Additional research work 7 demon- 
strated that the modal damping of composite structures depends also on the structural 
shape and structural deformation (mode shapes). This work also suggested that properly 
designed composite structures can provide significant passive damping, and they maj r fur- 
ther improve the dynamic performance and fatigue endurance by attenuating undesirable 
elasto-dynamic phenomena such as structural resonances, overshooting, and long settling 
times. The previous studies have also demonstrated that any increase in damping typically 
results in decreased stiffness and possibly strength, therefore, any tailoring of the compos- 
ite material for optimal damped response will be based on trade- offs between damping, 
stiffness, and strength. 

Reported work on the optimization of composite damping is mostly limited on local 
laminate properties instead of structural ones 8-9 . The authors have reported work on the 
optimization of the transient and forced dynamic response of composite structures 10 11 . 
This work is limited to a single objective function and laminate tailoring only, but il- 
lustrates that in order to realize full benefits from the damping capacity of composite 
materials, integrated methodologies for the optimal design of composite structures should 
be developed entailing: (1) multiple objectives to effectively represent the array of compet- 
ing design requirements; (2) capability for tailoring the basic composite materials and/or 
laminate; (3) capability for shape optimization; and (4) design criteria based on the global 
static and dynamic response of the composite structure. 

This paper presents the development of such a formal design method which allows 
the optimization of composite structures for optimal damped dynamic performance based 
on multiple objectives. The proposed design objectives are minimization of resonance 
amplitudes (or maximization of structural damping), minimization of structural weight, 
and minimization of material cost. The structural damped dynamic response is simulated 
with finite element analysis. A preliminary evaluation of the methodology on the struc- 
tural optimization of a cantilever composite beam is presented. The results quantify the 
importance of structural damping in improving the dynamic performance of composite 
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structures, and illustrate the suitability of multiple objectives in the optimal design of 
composite structures. 

Damped Structural Dynamic Response 

Assuming that finite element discretization has been applied, the dynamic response 
of a structure which is excited by a force P(t) is expressed by the following system of 
dynamic equations: 

[M]ii + [C]\x + [K]n = P(t) (1) 

where u is the discretized displacement vector. In the case of laminated composite struc- 
tures, the stiffness, damping, and mass matrices, [K], [C], and [M] respectively, are synthe- 
sized utilizing micromechanics, laminate, and structural mechanics theories representing 
the various material and structural scales in the composite structure. 

The related theories for this multi-level simulation of structural composite damping 
are described in refs. 1,2, and 7. Analogous theories are utilized for the synthesis of other 
mechanical properties 12 . At the micromechanics level, the on-axis damping capacities of 
the basic composite material systems are calculated based on the constituent properties, 
material microstructure, fiber volume ratio (FVR), temperature, and moisture. The off- 
axis damping capacities of the composite plies are calculated at the laminate level, and 
the local laminate damping matrices are predicted based on on-axis damping values, ply 
thicknesses, and laminate configuration. The damping contributions of the interlaminar 
damping layers due to interlaminar shear stress are also incorporated. 

The structural modal damping is synthesized by integrating the local laminate damp- 
ing contributions over the structural volume. The modal specific damping capacity (SDC) 
of the n-th vibration mode is: 


^ J A SW Ln dA 

j A w m Ja 


( 2 ) 


where: A is the structural area; SWi n and Wi n are respectively the dissipated and maxi- 
mum stored laminate strain energy distributions of the n-th mode per unit area per cycle. 


5 



Utilizing the finite element discretization scheme proposed in ref. 7, the modal SDC is 
finally related to the finite element damping and stiffness matrices, \C t ] and [A" e ] respec- 
tively: 


= 


nef 

y i u . ^ 

? . 2 Uein 

t=l 
TIC l 


’[Ceiju* 


(3) 


£ §u ein T [Ar e i]u«i n 

t=l 

where, nel is the total number of elements, and u e i n the nodal displacements of the i-th 
element corresponding to the n-th vibration mode. The dynamic response of the' structure 
is simulated based on modal superposition. The dynamic system in eq. (1) is transferred 
to the p x p modal space via the linear transformation u = [$]q. Assuming proportional 
damping and utilizing the first p modes, the following set of p uncoupled dynamic equations 
results: 


[m]q + [c]q + [*Jq = p(t) 


(4) 


The frequency response of the structure would be, 


«(/) = l*Mf) 


(5) 


where q(/) is the frequency response in the modal space induced by a harmonic force 
P3in(27r/f) of frequency /. A typical frequency response function (FRF) is shown in 
Fig. 1. The dynamic amplitude Uj at a frequency / is the complex magnitude of the 
displacement 

VM) = M/)| (6) 

The resonance amplitude Uj n at the n-th vibration mode is: 

v;„ = i«, •(/*.)! (7) 

where j^ n is the n-th damped natural frequency. The dynamic and resonance amplitudes 
in eqs. (6) and (7) are used as dynamic performance measures. The resonance amplitudes 
are primarily related to modal damping, and they decrease as the damping increases. 
Increased modal stiffness will also reduce the respective resonance amplitudes. 
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Multi-Objective Optimal Design 

As previously explained, damping is one of the factors that control the dynamic per- 
formance of composite structures near resonances. High damping values result in improved 
vibration control and fatigue endurance. However, increases in composite damping may 
typically result in stiffness/strength reduction and/or mass addition. Proposed design 
methods for optimizing the undamped dynamic performance of composite structures will 
fail to utilize the full potential of composite materials and may increase the dynamic am- 
plitudes near resonances. 

The minimization of weight and material cost are also included in the objectives. The 
material cost is a crucial factor, restricting in many cases the use of composite materi- 
als. The distinction between weight and material cost is also stressed. Fiber reinforced 
composites are non-uniform materials, hence, weight minimization does not correspond 
to material cost minimization. The design of composite structures for optimal dynamic 
performance is a multi-objective task, and may be best accomplished as the constrained 
minimization of multiple objective functions. 

A constrained multi-objective problem involving minimization of l objective functions 
is described in the following mathematical form: 

min{fi(z),F 2 (z) r .,F((z)} (8.1) 

subject to lower and upper bounds on the design vector z, inequality constraints G(z), 
and equality constraints H(z): 

z L < z < z u (9.1) 

G(z) < 0 (9.2) 

H(z) = 0 (9.3) 

In the rest of the paper, upper and lower values are represented by superscripts L and 
U respectively. Individual minimizations of each objective function subject to constraint 
set (9) will result in a set of l suboptima Fi = (Ff , F 2 , ..., F/), F 2 = (Fi , F £, ..., Fi), and 
so forth. These suboptima define a target point F* = (i^*, jF 2 *, ..., Ff) in the objective 
function subspace which usually lies in the infeasible domain. Hence, a candidate solution 
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may be obtained by finding a feasible point F = in the objective function 

subspace as closely as possible to the target point F*. This is equivalent to minimizing 
the distance between points F and F* subject to constraints (9), that is: 


minJIF; - Ff|| 


( 10 ) 


where the symbol j| || implies an admissible metric. In the present study, a weighted 
Euclidean metric is chosen for its simplicity and physical meaning. Other metrics may 
also be utilized, but in general, they are expected to result in different solutions. Hence, 
the original multi-objective problem becomes a constrained minimization problem of the 
following scaled objective function: 


min 


x> 

i=l 


{Fi-K’f 

F* 2 

» 


( 11 ) 


subject to constraints (9). The weighting coefficients are represented with u,-. 

In the pre sent method, the design objectives include minimization of: (1) the maxi- 
mum resonance amplitude (F\ = max{UJ n })] (2) the total structural weight W ( F2 = W)\ 
and (3) the material cost represented by the average cost of fibers per unit area Pf 
(i *3 = Pf). The resonance amplitudes are related to modal stiffness, modal damping, 
and modal mass, hence, the objective function F\ will tune the modal stiffness, damping 
and mass. Alternatively, F\ may represent modal damping values. The explicit maxi- 
mization of modal damping may be preferred in the case if an a-priori unknown dynamic 
excitation. The use of fiber cost as a measure of the total material cost is justified in 
view of the very high cost of fibers compared to the cost of matrix. The design vector 
includes fiber volume ratios (FVRs), ply angles, and shape parameters. The proposed 
design criteria are then formulated in the following form: 


min{F u F 2 ,F 3 } 


( 12 ) 


subject to constraints (9.1), constraints on static deflections u B , 


u 8 < u 


u 


(13.1) 
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dynamic amplitudes (including dynamic resonance amplitudes), 


U 4f) < Uj (13.2) 

u; < (i3.3) 

and natural frequencies f n . 

f : L < f„ < f u (13.4) 

This constrained optimization problem is solved with non-linear programing. In the 
present paper the modified feasible directions non-linear programing method 13 is utilized. 

Application and Results 

The proposed design method was applied on the optimal tailoring of a cantilever 
graphite/epoxy composite beam. The initial beam shape (6in long, lin wide, and 0.2in 
thick) is shown in Fig. 2a. Typical properties of the composite material in room con- 
ditions are shown in Table 1. The SDCs of the fibers and matrix were backcalculated 
from the respective laminate SDCs reported in ref. 4. The beam was assumed to operate 
in room conditions. The composite beam is symmetric consisting of angle-plied sublam- 
inates 1, 2, and 3 in each side. All sublaminates have plies of equal thickness (O.Olin). 
Sublaminate 3 is at the center of the beam. Fig. 2b shows the laminate configuration 
((±#i) 2 /(±# 2 ) 2 /(± 03 ))j for the initial uniform beam thickness. The ply angles 6{ and 
FVRs kfi of each sublaminate, and the thicknesses hj at 0%, 30%, 60%, and 100% (tip) 
of the beam span were the design variables. The thickness at other sections of the beam 
was interpolated using a cubic spline fit. It was further assumed that when the thickness 
of a cross-section was reduced, the inner plies of Sublaminate 3 would first drop, and if 
necessary, the inner plies of Sublaminates 2 and 1 would drop subsequently. A central 
ply was dropped when its thickness was reduced to less than O.OOSin. In such case the 
thickness of the adjacent ply was increased by an equal amount to avoid discontinuity. 
In the opposite case of increased thickness, inner plies of O.Olin thickness were added to 
Sublaminate 1. 
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A finite element mesh of 55 nodes and 80 specialty triangular plate elements was 
utilized. A static transverse out-of- plane (y-axis) force of 5 lbs/in, and a transverse out-of- 
plane harmonic force of 0.1 lbs/in amplitude were applied at the tip of the beam. Under 
this type of dynamic loading, the maximum resonance amplitude at the tip corresponds to 
the first mode (first out-of-plane bending) (Fi = U ^). Because a single composite system 
w r as assumed, the material cost was replaced by the average fiber volume per unit area Vf, 
( Fs = Vf). The weighting coefficients in eq. (11) were set equal to unity. Constraints were 
imposed on the transverse static deflections and on the transverse resonance amplitudes 
of the first five modes at the tip of the beam. Additional lower bounds were imposed on 
the first two natural frequencies. 

u iy < 35 10 _3 m (14.1) 

Uy n < 0.45m (14-2) 

fi > 400 Hz and / 2 > 1200Hz (14.3) 

The following upper and lower bounds were imposed on the design variables: 

-90.0° < 6i < 90.0° (15.1) 

0.01 < k fi < 0.70 (15.2) 

0.04m <hj < 0.4m (15.3) 

Table 2 shows the initial reference design, the three suboptimal designs with each 
objective function individually optimized, and the resultant optimal design of the multi- 
objective optimization. The unidirectional beam was selected as the initial reference design 
because it provides the higher bending rigidity. The predicted maximum static deflection, 
as well as the resonance amplitudes at the tip, natural frequencies, and modal SDCs of the 
first 3 modes are also shown. As seen in table 2, the initial design violates both frequency 
constraints, but all optimized designs have natural frequencies in the feasible frequency 
domain. All optimized designs have non-uniform thickness, being thick at the proximal 
end and thin at the distal end. Side views of the optimized thickness shapes for each 
optimization case are plotted in Fig. 3. Also, the static deflection constraints were active 
only in the minimum- weight design. 
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The optimized designs differ in many aspects. The minimization of the resonance tip 
amplitude ( minFi ) decreased the resonance peak by 65% with respect to the unidirec- 
tional beam. This reduction was attained by a combination of high modal stiffness and 
damping. The shape optimization resulted in increased weight and material cost (12%) 
and redistribution of the material. The near ^27° ply-angles in the outer two sublami- 
nates contribute to high flexural damping and increased in-plane shear stiffness, while the 
high FVRs provide out-of-plane flexural stiffness. The inner ply primarily contributes to 
damping. 

The minimum weight design decreased the weight of the beam by 23%, the 

material cost by 31%, but also increased the resonance amplitude by 12% with respect to 
the initial design. The weight reduction was primarily accomplished by material redistri- 
bution, lower ply angles in the outer sublaminate, and lower FVRs in all sublaminates (the 
fibers have higher density than the matrix). In fact, the inner sublaminate was virtually 
reduced to pure matrix resulting in a sandwich type laminate configuration. The first 
modal damping has moderately increased. 

The minimum cost design ( minFz ) decreased the material cost by 50%, the structural 
weight by 15%, and increased the first resonance amplitude by 11%. Interestingly, the 
ply angles of the two outer sublaminates are ±3.2° and ±3.7° respectively, an almost 
unidirectional configuration, and the FVR’s were reduced. The inner sublaminate was 
again reduced to pure matrix. This laminate tailoring provides the ultimate utilization of 
the reinforcing fibers. Both natural frequency constraints are active, and the first modal 
damping of the optimized beam is low. 

Clearly, all individual minimizations of each objective function failed to produce simul- 
taneous improvements in all objectives. In contrast, the resultant multi-objective optimal 
design seems to blend many of the merits of the individual single-objective optimizations, 
since all objective functions were significantly reduced compared to the initial design. 
Specifically, the first resonance amplitude was decreased by 60%, the structural weight by 
10%, and the material cost by 50% . This clearly illustrates the superiority and suitability 
of the multi-objective optimal design for composite structures. The optimum laminate 
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configuration is composed by an outer ^26° angle-ply sublaminate of high FVR providing 
all the stiffness, and two inner sublaminates of virtually pure matrix. It seems likely, that 
the reduction of Sublaminate 2 to pure matrix was driven by the requirements of high 
damping and low material cost. Indeed, the present optimal design exhibits the higher 
value for the first modal SDC, and material cost reduction equal to that obtained by the 
minimum cost design. This suggests the possibility of utilizing hybrid laminates, and layers 
of high-damping viscoelastic materials for further improvements. 

The frequency response functions at the mid-point of the free-edge of the initial and 
optimized beams are shown in Fig. 4. The multi- objective optimum design produced 
the best FRF, while the minimum cost and minimum weight designs the most inferior. 
This suggests that the incorporation of composite damping was crucial in obtaining these 
significant improvements in all objective functions illustrating, in this manner, the signifi- 
cance of composite damping in the design of high- performance, light-weight, and low-cost 
composite structures. Optimal design methodologies neglecting the damping capacity of 
composite materials and its controllable anisotropy may lead to structures with inferior 
dynamic performance near the resonance regimes, hence, they may be unsuitable for opti- 
mizing the dynamic performance of composite structures. 

Summary 

An optimal design methodology for optimizing the dynamic performance of composite 
structures based on multiple objectives was proposed. The effects of composite damping on 
the dynamics of composite structures were incorporated. The design objectives included 
minimization of resonance amplitudes (or maximization of modal damping), minimization 
of structural weight, and minimization of material cost. Performance constraints were 
imposed on static displacements, dynamic resonance amplitudes, and natural frequencies. 
The described method has been integrated into an in-house research code (STAT) [14]. 

The method was applied on the optimal design of a cantilever composite beam. The 
multi-objective optimization was proved superior in optimizing the competing require- 
ment for high damping, low weight, and low material cost, and resulted in significant 
simultaneous improvements in all objective functions. Optimizations with single objective 
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functions did not improve all objectives. The obtained results for this primitive structure 
also illustrated that the damping capacity of composite materials is an important factor in 
designing light-weight, low-cost composite structures of improved dynamic performance. It 
was further shown that proposed methods for optimizing the undamped dynamic response 
of composite structures may fail to improve, or may deteriorate, the dynamic performance 
near resonance, as they neglect the tailorable damping capacity of composite materials. 

Overall, the application of the method on this primitive composite structure appeared 
very encouraging. Additional applications of the method on more complex structural 
configurations will be performed in the near future. 
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Table 1. Mechanical properties of HM-S/epoxy system. 
(70°F and 0% moisture) 


Epoxy HM-S Graphite 


E m — 0.500 Mpsi 
(3.45 GPa) 
G m = 0.185 Mpsi 
(1.27 GPa) 
Vw = 10.30 % 

= 11.75 % 

Pm =0.0440 lb/in s 


Efu — 55.0 Mpsi 
(379.3 GPa) 
Ef 22 = 0.9 Mpsi 
( 6.2 GPa) 
G fi2 — 1.1 Mpsi 
( 7.6 GPa) 
v fl2 = 0.20 
ipfu = 0.4 % 

V>/22 = 0.4 % 

V>/12 = 0.4 % 
p m =0.0703 lb /in 3 
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Table 2. Initial and optimum designs. 


Initial Single-Objective Designs Multi-Objective 

Design (min F\) (min F 2 ) (min Fs) Design 


Thickness, (in) 


0% span 

0.2 

0.2529 

0.2090 

0.2022 

0.2584 

30% span 

0.2 

0.2313 

0.2046 

0.2014 

0.2343 

60% span 

0.2 

0.2572 

0.1573 

0.2011 

0.2329 

100% span (tip) 0.2 

0.1053 

0.0400 

0.1426 

0.0400 

Ply Angles, 

(degrees) 





e 1 

0.0 

-26.28 

-13.87 

3.213 

-26.46 

e 2 

0.0 

-27.80 

-45.34 

3.678 

-39.89 

h 

0.0 

19.99 

23.17 

11.80 

-58.75 

Fiber volume ratios 





kf i 

0.50 

0.632 

0.530 

0.482 

0.605 

k f 2 

0.50 

0.700 

0.478 

0.136 

0.010 

kf 3 

0.50 

0.116 

0.010 

0.010 

0.010 

Objective Functions 





Fi (in) 

0.3977 

0.1361 

0.4467 

0.4399 

0.1580 

F 2 (lbs) 

0.0677 

0.0758 

0.0520 

0.0576 

0.0609 

Fs(in 3 /in 2 ) 

0.0990 

0.1114 

0.0685 

0.0492 

0.0486 

Max. y-axis 

Static deflection (tip), 

(10- 3 in) 



u ty 

19.70 

21.24 

35.00 

24.38 

30.69 

Natural Frequencies, (Hz) 





h 

387.4 

400.2 

485.6 

400.0 

421.2 

h 

1095.3 

1462.3 

1971.0 

1200.0 

1362.2 

h 

1549.7 

2321.7 

2010.4 

1367.0 

2177.1 

Resonance y 

-axis Amplitudes (tip), 

(in) 



TJ T 

u vl 

0.3977 

0.136 

0.450 

0.443 

0.1577 

ul 

0.0000 

0.000 

0.002 

0.000 

0.0000 

ul , 

0.0000 

0.005 

0.071 

0.000 

0.0102 

Modal SDC’ 

s, (%) 





0i 

0.604 

1.881 

0.807 

0.654 

2.238 

02 

5.592 

3.374 

1.387 

5.251 

3.769 

03 

3.503 

2.287 

0.869 

2.992 

2.533 
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AMPLITUDE 



FREQUENCY <M LAMINATE CONFIGURATION. 

FIGURE 1. - TYPICAL FREQUENCY RESPONSE FUNCTION. FIGURE 2. - CANDIDATE COMPOSITE BEAM. 



FIGURE 3. - INITIAL AND OPTIMUM THICKNESS DISTRIBUTIONS. FIGURE A. - FREQUENCY RESPONSE FUNCTIONS OF THE INITIAL 

AND OPTIMAL DESIGNS. 
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